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Abstract 

1^ • We present numerical visualizations of Ricci Flow of surfaces and 3-diniensional manifolds of 

On . revolution. Riccijrot is an educational tool which visualizes surfaces of revolution moving under 

I I Ricci flow. That these surfaces tend to remain embedded in M'^ is what makes direct visualization 

VO I possible. The numerical lessons gained in developing this tool may be applicable to numerical 

^-r ' simulation of Ricci flow of other surfaces. Similarly for simple 3-dimensional manifolds like the 

^3 ' 3-sphere, with a metric which is invariant under the action of S0{3) with 2-sphere orbits, the 

^ . metric can be represented by a 2-sphere of revolution, where the distance to the axis of revolution 

j^ ! represents the radius of a 2-sphere orbit. Hence we can also visualize the behaviour of such a 

H I metric under Ricci flow. We discuss briefly why surfaces and 3-manifolds of revolution remain 

IJ ' embedded in M^ and R'* respectively, under Ricci flow and finally indulge in some speculation 

about the idea of Ricci flow in the larger space of positive definite and indefinite metrics. 



X 



1 Introduction 

Understanding Ricci flow has become fundamental to attempts to prove Thurston's ge- 
ometrization conjecture (|3], |H], [H], [IHl)- Since Ricci flow is geometric by nature, it 
is natural to ask whether it can be visualized. Unlike mean-curvature and related flows, 
which act in the flrst instance as a force pushing on a surface and therefore tend to keep 
surfaces embedded in Euclidean space, Ricci flow acts directly on the metric of the surface, 
tending not to preserve embeddedness. A number of interesting results have been obtained 
by restricting to classes of metrics of revolution, since such symmetries are preserved under 
Ricci flow and the metric depends on considerably fewer parameters in such cases f|21j. 
j2j, Pn])- We will discuss for this class of metrics, why Ricci flow does preserve embed- 
dedness in Euclidean space (PJ). Moreover, the key phenomenon of neck pinching occurs 



naturally for metrics of revolution. Good estimates are available to understand the limit- 
ing behaviour of metrics as they become singular under neck pinching in (|2I], EI)- Our 
visualization approach enables us to generate interesting pictures and information about 
neck pinching. In ^2J, it is shown that any metric on the 2-sphere converges to the round 
metric under Ricci flow, so neck pinching occurs only in higher dimensions. 

Note that there are many interesting questions about Ricci flow for 3-dimensional man- 
ifolds. There are a limited number of examples where one knows completely how the 
metric evolves ([5|)- It would be very useful to have further examples where Ricci flow 
continues for infinite time, without developing singularities. For example, there are im- 
portant constructions (Dehn surgery) of metrics with strictly negative sectional curvature 
on 3- manifolds ([H], |H1)- But even in this case, it is not clear whether singularities will 
form under Ricci flow. In fact, no examples are known where a metric with all sectional 
curvatures negative flows to a metric with positive sectional curvature at some points. 

2 Preserving the Embedding into Euclidean Space 

In this section, we give a quick discussion of why n-dimensional spheres of revolution 
{K2 = 1 - see Section 15. Ij) . which are initially isometrically embedded into Euclidean 
(n + l)-dimensional space, remain isometrically embedded, so long as Ricci flow has a 
smooth solution. 

2.1 Characterising when Isometric Embeddings exist 

We observe a simple extension of [Hj. Note that in later sections, we will use extrinsic 
coordinates, which are useful to see how points of the surface or manifold move under 
Ricci flow when visualized in Euclidean space. Here it is more convenient to use intrinsic 
coordinates. So on S"", a metric is chosen of the form ds^ + a^{s)d6'^, where s is arc length 
along a geodesic joining the North Pole to the South Pole and 6 represents the coordinates 
in a hypersphere of radius q;(s), which is the orbit of an isometric action of SO{n) on S"', 
where the orbit space is an arc. We assume that s varies between and L. 

In 0, it is shown that a necessary and sufficient condition that S"^, with a metric of this 
form, can be isometrically embedded in M^ as a surface of revolution, is that the integral 
of Gaussian curvature over any "polar cap" is positive, where such a region consists of all 
points where either 0<s</!;orA;<s<L, for some k with < k < L. Note also that 
a surface of revolution in M^ is invariant under an isometric action of 5*0 (2) by rotation 
about the x^ axis. Similarly we define a manifold of revolution in M^"+^^ as being invariant 
under the isometric action of SO{n) by rotation about the x^ axis. Finally we follow |9|, by 
requiring that any manifold of revolution ]R*^"+^) has at most one component of intersection 
with any hyperplane of the form x^ = k for any constant k. Equivalently, the manifold 
is obtained by rotating the graph of a function of one variable around the x^ axis, by the 
action of SO{n). Note that for visualization purposes, it is not essential but is certainly 
convenient to have this restriction. 



Theorem 1 Suppose that a metric of the form ds'^ + a'^{s)d9'^ is chosen on S^. Then there 
is an isometric embedding as a manifold of revolution in M.^"^^^^ if and only if the same 
metric, viewed as on S'^ , can be isometrically embedded into M^ as a surface of revolution. 

The proof of this theorem is very easy and is left to the reader. 

2.2 Ricci Flow preserves Embeddability of Manifolds of Revolu- 
tion 

An observation about the characterisation in [^ is as follows. For a polar cap, Gauss 
Bonnet shows that the sum of the integral of the Gaussian curvature of the cap and the 
integral of the curvature of the boundary curve C is it. Therefore positivity of the first 
integral is equivalent to the second integral being strictly less than tt. Note that C is the 
orbit of SO (2) acting on a plane parallel to the x'^x^ plane and so is a standard round 
circle in this plane, with centre on the intersection of the plane with the x^ axis. Hence the 
integral of its curvature in M^ is exactly vr and the direction of its curvature at each point 
points towards the x^ axis. But then to compute the curvature in the surface of revolution, 
one has to project onto the tangent space of this surface. Hence the projected curvature 
vector will be shorter, unless the tangent plane is vertical at every point of C. So this 
shows that the condition on the polar caps remains true as Ricci flow proceeds, unless at 
some time and for some value of s, with < s < L, |a'(s)| = 1. Therefore we need to 
show that |a'(s)| < 1 remains true, away from the poles, so long as Ricci flow produces a 
smooth solution ( , [01 , [IBl ) , to prove the following result. 

Theorem 2 // there is an isometric embedding of S*" as a manifold of revolution in R*^""^^) 
then the manifold remains isometrically embedded so long as Ricci flow gives a smooth 
solution. 

Proof Let v = a', where v and a are viewed as functions of s and time t, and as 
above, the derivative of a is taken with respect to s. By equation (16) of ,2j, t; satisfies 
the following evolution equation. 

- n — 2 n — !,_, n. 
OtV = Vss H VVs H ^(1 - V )v. 

Here we are using v^ to denote partial differentiation in the s variable. 

Now initially, \v\ = \as\ < 1 for all values < s < L, by assumption. Let t be a first 
value of time for which |f(s,t)| = 1, at some s not corresponding to the North or South 
pole of the sphere. If v{s, t) = 1, then this is a maximum value of a whereas if v{s, t) = —1, 
then this is a minimum value. 

So we see by the maximum principle, (see eg |H|) that this gives a contradiction, since 
1 — v^ = Vs = ioT both cases. This shows that the sphere does remain embedded, so long 
as it is smooth. 



3 Computational Formulation 



A two-dimensional surface of revolution of genus zero embedded in M.^ can be defined in a 
polar representation with coordinates x^ = p and x^ = 9 & [0, 27r[ by a metric of the form 
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where \Jm{p) has the direct physical interpretation as the radius from the axis of rotation. 
For a closed surface, we require 

m{ppoie) = 0. 

We have chosen PNorth_Pole — ^^^ PSouth_Pole — ^- When h{p) is a constant, \/h x p 
is the distance of a point with coordinates (p, 9) from the North Pole along a meridian. 
Smoothness at the poles demands 
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Note that the necessary and sufficient condition for isometric embeddability as a classical 
surface of revolution in E^ for iS^-invariant metrics on S"^ - the nonnegativity of the integrals 
of the curvature over all disks centred at a pole - is given in jU]. The generating curve 
(cross-section) of the surface of revolution is given by 
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with p e [0, n]. 

The Ricci tensor is 
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where primes denote partial differentiation with respect to p. 
Ricci flow satisfies 



at 



-2R 



■flU- 



(2) 



\pe]0,7r[) 



(3) 



:pG{0,7r}), 



The possible initial metrics Ricci_rot provides are given by 











sin p + C3 sin 3p + C5 sin 5p 

1 + 3C3 + 5C5 



(4) 



These surfaces are symmetric under reflections about a plane normal to their axis of rota- 
tion which passes through their centre of mass. This extra symmetry has proven useful in 
making the simulation faster, since it then suffices to compute on only half of the surface. 
Attempting to integrate using Equation Q and central finite differences in a purely 
explicit scheme leads, not surprisingly to numerical instability, particularly (but not only) 
at the poles. This instability can be reduced significantly by noting that 
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3.1 Filtering and Reparametrization 

The numerical instabilities are however persistent enough to require more potent measures. 
We did not try any implicit scheme, which may have alleviated the numerical instabilities, 
but would still have required very small time steps near the poles. 

Inspired by spectral methods |1], we first introduced a filter which consists of trans- 
forming to Fourier space (DFT), dropping shorter wavelength terms, and then transforming 
back. This process is facilitated by noting that the rotational and refiection symmetries of 
the family of surfaces (0)) are preserved by Ricci fiow. We can therefore write 

h{p) = 2^ hi cos 2ip 



i=0 



and 
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j=0 

where equality of course only applies to Nh = N^ = C)0. We actually choose N^ and A^^ 
to be much less than the number of points used in the finite difference scheme, and find 
that there is no need to implement the FFT. We need to ensure that (0) is satisfied at the 
poles, to be sure that the generating curve can be found. We do this by multiplying 
the computed values of y/m by 
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where K is a small positive constant, but large enough that this correction factor quickly 

becomes unity away from the poles. Setting K = has the unwanted side-effect that local 

numerical errors near the poles can lead to global changes in shape which do not satisfy 

the Ricci flow. 

By its very nature, the Ricci flow forces some parts of a surface to contract while others 

expand. This means that whatever set of nodes we have chosen for our flnite-difference 

scheme will effectively become less and less evenly distributed, creating further sources of 

numerical instability where they become too dense. The solution to this problem is to 

recall that any given surface corresponds to many different parametrizations, each with a 

corresponding metric. From a numerical point of view, the most pleasant parametrizations 

let h{p) be a constant (this corresponds to p being directly proportional to distance from 

a pole). Let 

p 



s=0 



be the distance from the North Pole. We have the freedom to reparametrize: 
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where p G [0, tt] both before and after. It is this reparametrization which, performed after 
every few iterations of the explicit integration and flltering steps, has made Ricci_rot 
stable enough to use. 

4 Use of the Software 

First, a quick introduction: Press "n" and then drag the mouse to choose an initial surface. 
Then hold down the up-arrow key to watch the flow until it stops. Once one has begun 
flowing, dragging the mouse will rotate the surface. 



4.1 Details 

Ricci_rot is an ANSI C program ^H] which uses the OpenGL standard ^\ to display 
graphics. It should therefore be portable, although it has only been tested on a Mac OS 
X platform. 

Upon launch, a window is opened with the image of a sphere with blue meridians and 
parallels. This initial shape can be deformed by holding down the left mouse button and 
dragging in any direction. The possible initial shapes belong to the family Q, where c^ is 
varied by horizontal, and C5 by vertical motion. The program will not accept all possible 
values of C3 and C5, so, dragging very far in any direction, one may flnd that the shape 
ceases changing. 



To examine the shape, pressing "f" (for "flow") will put the program into flow mode. 
The meridians and parallels change colour. Now, dragging the mouse rotates the surface. 
To make any more changes, pressing "n" (for "new shape") will put the program back in 
its initial mode. The meridians and parallels change colour back to blue, and dragging the 
mouse changes the shape. 

When one has chosen a shape one wishes to flow, it is enough to begin pressing the 
up-arrow key. Usually this key will repeat if held down. Then one can see the surface flow 
continuously. Since the program is now in flow mode, dragging the mouse will rotate the 
surface. 

Once the program has met with a numerical instability, the flow will stop, and the 
parallels and meridians will appear black. Any ripples which may appear on the surface 
at this stage are a result of the numerical instability. 

At any time, pressing the left-arrow or right-arrow keys will rotate the surface anti- 
clockwise or clockwise respectively. 

The down-arrow key flows the surface backwards in time. This evolution is highly 
unstable. One should not expect to be able to flow for very long. 

Pressing "m" at any stage will change the display mode from showing the surface to 
showing the components of the metric {gn = h in green and (722 = ''^ in blue) and the 
surface's cross-section in white. To revert to displaying the surface in E^, press "s". 

There is a bright spotlight. This can be turned on or off by pressing the right mouse 
button. A menu will appear, from which one can choose to turn the spotlight on or off 
(or change to and from "flow" and "new shape" modes). It can be rotated (in a rather 
non-intuitive way) by dragging the mouse with the middle mouse button held down. 

4.2 Examples 

The first example is of a dumbbell shape, given by Equation Q with C3 = 0.766 and 
C5 = —0.091. This surface's flow is illustrated in Figure ^ with time step dt = 0.01. The 
flow of its cross-section is provided in Figure |2l time step dt = 0.002. It is interesting to 
compare this with Figure 4 of jnj, which is an illustration of the mean-curvature flow of 
the cross-section of what is initially a dumbbell shape. 

Figures 01 IH El and IHl were all generated using the initial surface given by C3 = 0.021 
and C5 = 0.598, and with the time step dt = 0.002. The values oi gn = h {h{p) becomes a 
constant after reparametrization) are 1, 1.027938, 1.000526, 0.936608, 0.843907, 0.729106, 
0.601463, 0.475008, 0.365054, 0.278860, 0.213761 and 0.163754. 

The behaviour is in both cases what one expects - flow towards spheres of constant 
positive Gaussian curvature [7j. See also Chapter 5 of jH], in which it is shown that any 
solution of the unnormalized Ricci flow on a topological S"^ shrinks to a round point in 
flnite time, following [.7|- Note that the "modifled" Ricci flow in jl6j, for the special case of 
surfaces of revolution embedded in Euclidean space, could be visualized in a similar way. 
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Figure 1: Deformation of a dumbbell-like surface of revolution under Ricci flow. The pictures 
have been taken at equal time intervals and are drawn to the same scale. 




Figure 2: Deformation of tire cross-section of a dumbbell surface of revolution under Ricci flow. 
The axis of rotation is horizontal. The pictures have been taken at equal time intervals. The 
initial and final curves are both shaded black. 

5 Three-Dimensional Manifolds of Revolution 



5.1 A Line cross a Surface of Revolution of Constant Curvature 

Let K2 be the Gaussian curvature of a general abstract Riemannian surface of revolution, 
where K2 is any real constant. 
The 3-nianifold metric is 
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where p = x^ plays the role of a latitude and 6 = x'^ the role of a latitude on the abstract 
Riemannian surface of revolution. 

The non-zero elements of the Ricci tensor are 
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and therefore unnormalized Ricci flow satisfles 
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Figure 3: Deformation of a surface of revolution under Ricci flow. The pictures have been taken 
at equal time intervals and are drawn to the same scale. The cross-sections are plotted in Figure 
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Figure 4: Deformation of a surface of revolution under Ricci flow. The pictures have been taken 
at equal time intervals and are drawn to the same scale. These pictures are the continuation of 
the evolution of the surface begun in Figure |31 The cross-sections are plotted in Figure [3 
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Figure 5: Deformation of the cross-section of a surface of revolution under Ricci flow. The axis 
of rotation is horizontal. The pictures have been taken at equal time intervals. The initial and 
final curves are both shaded black. 
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The scalar curvature is 
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the sectional curvatures are 
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5.2 An Example 
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We will study the pinching behaviour under unnormalized Ricci flow of the initial 3- 
manifold with metric given by 
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h{p) = 1 and K2 = 1 
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at t = 0. We denote the time at which pinching occurs (m(0) = 0) by t = T. 
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Figure 6: Evolution of the metric component 522 = rn{p) of a sm'face of revolution under the 
combined action of Ricci flow and reparametrization - meaning that the corresponding curves 
gii = h are independent of p. The curves have been plotted at equal time intervals. The initial 
and final curves are both shaded black. 

5.2.1 Qualitative Behaviour 

The most direct approach to simulating the flow determined by Equations IHl and ^1 is to 
use an explicit finite-difference formulation. Doing this, we found that the system is (not 
surprisingly) highly unstable. We identified two modes of instability. One was of a short 
wavelength - of the type k x (—1)* where i is the spatial grid index and k some constant. 
This type of instability could be removed by simply performing the simulation in Maple 112] , 
using several hundred digits' precision. The other type of instability showed itself as a long- 
wavelength smooth perturbation which diverged as the number of time-steps was increased. 
We were not able to eliminate this second source of instability. Practically, it meant that 
we were restricted to fairly large time-steps. To be more precise, we stepped between 
t = 0.000025, 0.000050, 0.0000625, 0.00006875, 0.000071875, 0.000075000, 0.00007578125, 
0.00007656250, 0.000076953125 (shown only in Figure El) and 0.000077343750. Since we 
are therefore not able to claim that the method converges to the exact solution, we have 
chosen to speak of the data from this simulation as being only qualitatively correct. 

It is nonetheless informative to consider Figures [7| |H1 and IHl 

Figures |H1 and El would suggest that pinching occurs before t = 0.000077343750. At 
t = 0.000077343750, values of m at nodes near p = are negative. Since Equations El 
no longer make sense for negative m, an interval around p = has essentially been cut 
out. This would of course correspond to the view that pinching creates two independent 
geometric bodies, as suggested in Figure El Our single coordinate p should then be replaced 
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Figure 7: Qualitatively correct evolution of the metric component gn = h{p) of a 3-manifold of 
revolution under Ricci flow. The curves have been plotted at unequal time intervals. 

by two new coordinates, one for each body, and values of p which can no longer be attributed 
to the one or the other body should then be discarded as being "non-geometric" . It is indeed 
interesting to consider what might happen after pinching, retaining the coordinate p and 
dropping the condition that the metric should be positive definite. It is however important 
to note that the negative values of m seen in our qualitative numerical results cannot be 
trusted, since they are an artifact of Euler's method being applied with a step-size so large 
that a singularity was passed over. What we do claim is that the cross-sections plotted 
in Figure El do still have some meaning, if one does not attempt to read anything into the 
missing points at the caps of the two bodies which resulted from pinching. 
We return to this question in Section IHl from another point of view. 

5.2.2 Asymptotic Scaling 

A second approach to studying pinching involves series expansions of the metric compo- 
nents h and m. We have expanded them to tenth order in p (even terms only) and first 
order in t: 
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h{p,t) = y^{hi + hit) p'' and rn{p, t) = / ^(m^ + mj t) p\ 
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Substituting into Equations Q and (fTUI) we get, by equating coefficients (and assuming 
higher-order terms are zero) equations for the quantities hi and mi. These equations are 
quite cumbersome. We have used Maple jI2] both to derive the equations and also to 
generate C code jT3] to allow us to compute the flow as rapidly as possible. The C code 
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Figure 8: Qualitatively correct evolution of the metric component 922 = Tn{p) of a 3-nianifold 
of revolution under Ricci flow. The curves have been plotted at unequal time intervals. 

was modified to utilize the long double numerical data type, and the code was run on a 
DEC alpha (the name of a RISC microprocessor) machine with true 64 bit floating-point 
arithmetic (the long double type). We used Euler's method with a very primitive form 
of adaptive step-size estimation, flnding these to be satisfactory for our purposes. 

Our data are consistent with the qualitative results of the previous section. 

We are able to use these series expansions to make some conjectures concerning the 
behaviour of various quantities as t approaches T. The following scaling laws are purely 
empirical. They are the result of curve-fltting to the data we have, in particular for the 
seven points, at t = 0.000079300, 0.000079310, 0.000079320, 0.000079330, 0.000079340, 
0.000079345 and 0.000079350. The step-sizes used were signiflcantly shorter than the 
differences between the times of these representative points. Figures HU] and [TT] illustrate 
this fltting process. It is clear from these illustrations that what we believe to be asymptotic 
behaviour sets in fairly late. One may of course ask if we have fltted to data which is late 
enough in the flow to truly warrant being called asymptotic. At this stage, we can only 
present the data we have. 

All of the following apply to p = 0: 
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1.409 X (0.0000793514 -t) ^-^^^ 
1.705 X (0.0000793529 - t)-o-235 
0.570 X (0.0000793515 - t)-^-°^^ 
-1.142 X (0.0000793513 - t)-o-826 
0.698 X (0.0000793514 - t)''^-^^^ 
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It would therefore appear that T ^ 0.0000793514. 
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Figure 9: Qualitatively correct deformation of the cross-section of a 3-manifold of revolution 
under Ricci flow. The axis of rotation is horizontal. The pictures have been taken at unequal 
time intervals. 



Note that this curve- fitting problem is ill-conditioned, although Figure [TT] does suggest 
that our fits are as accurate as one could hope for. However, h appears to diverge at a 
significantly later time than our estimate for T, and we do not find that our curves of best 
fit for m and K{h., c) obey m x K{h^ c) = K2 exactly at p = 0. 

We do find that the second derivative of the cross-section y{x) (see Equations Ej) at 
X = diverges as t — *> T, indicating that the neck becomes sharper rather than longer as 
the surface approaches pinching. This makes sense, since at t = T we expect to have two 
abutting caps, as illustrated in Figure IHl 

These results need to be reproduced independently (perhaps via the DeTurck flow [TU] ) 
to ascertain their validity. We present them on an "as is" basis, freely admitting that 
we are not even in a position to give an error analysis, but in the hope that even such 
preliminary data may inspire or support some new analytical attack. 

It is perhaps appropriate to quote Kenko at this point (from essay 82 of ^14j): 

Leaving something incomplete makes it interesting, and gives one the feeling 
that there is room for growth. Someone once told me, "Even when building 
the imperial palace, they always leave one place unfinished." 



One is tempted to think of level set methods (see [^ for a general introduction) as a 
possible computational alternative to what we have tried. Level set methods have been 
successfully applied to mean-curvature flow [HI, and indeed their strength lies in their ability 
to cope naturally with changes in topology of an evolving surface, as pinching involves. 
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Figure 10: The evolution of h with time for a 3-manifold of revolution. The points are numerical 
data. The curve is a fit to the asymptotic behaviour of h near pinching. 

We feel that some reformulation of the problem, perhaps along the lines of what fol- 
lows, may facilitate a level set formulation of Ricci flow, and that this is without doubt a 
worthwhile focus for further work. 



6 Speculation 

The system of differential equations and ()1()|1 above makes sense when the metric tensor 
is not necessarily positive definite, so long as smoothness occurs at points where m = 0. So 
we are lead to speculate about the idea of Ricci flow in the larger space of positive definite 
and indefinite metrics. 

The advantage of such an approach might be that instead of the manifold pinching, 
one might allow that part of the manifold has a positive definite metric and part has an 
indefinite metric, after a singularity occurs. In this case, the computational problem might 
become more tractible. The reason is that the flowing family of indefinite metrics would 
always have a fixed domain of definition. 

In particular, at a neck pinch, the initially Riemannian metric for a rotationally sym- 
metric sphere should "jump" to a metric with a positive definite part on two polar caps 
and an indefinite part on an equatorial band. The indefinite part will expand and the 
positive definite parts contract until the metric on the whole space becomes indefinite. 
Further numerical experimentation might give a better idea of a suitable formulation for 
this process. 

Note that Perelman ^Hl has indicated how his notion of entropy can be used as a 
potential function for Ricci flow. Extending this approach to the larger domain of positive 
definite and indefinite metrics may lead to a true level set formulation of Ricci flow, in this 
space of metrics on a fixed manifold. Note that at singularities where some components of 
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0.0000793 



0.00007932 



0.00007934 



Figure 11: The evolution of h with time for a 3-manifold of revokition close to pinching. The 
points are numerical data. The dark curve is the fit 1.705 x (0.0000793529- t)"°-^^^. The 
lighter curves correspond to other attempts at curve-fitting, in particular in an attempt to match 
the pinching time of T = 0.0000793514. These are 1.294 x (0.0000793514- 1)-°-24, 1.586 x 
(0.0000793514-t)-0-23, and 5.389 x (0.0000793514-t)-0-i^ 

the Ricci tensor become infinite, one would expect the level sets to fatten ( see e.g. P). 
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